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The anomalous proximity effect between a d-wave super- 
conductor and a surface layer with small electronic mean free 
path is studied theoretically in the framework of the Eilen- 
berger equations. The angular and spatial structure of the 
pair potential and the quasiclassical propagators in the inter- 
face region is calculated selfconsistently. The variation of the 
spatially-resolved quasiparticle density of states from the bulk 
to the surface is studied. It is shown that the isotropic gap less 
superconducting state is induced in the disordered layer. 

PACS numbers: 74.50.+r, 74.80. F, 74.72.-h 



I. INTRODUCTION 

There is continuing experimental evidence that the be- 
havior of high temperature superconductors (HTS) can 
be understood in terms of the d-wave pairing scenario, 
rather than in the conventional s-wave picture. On the 
other hand it is well known that the d-wave order param- 
eter is strongly reduced by electron scattering at impuri- 
ties and therefore can be formed only in clean materials. 
However, the condition of clean limit is not fulfilled in the 
vicinity of the grain boundaries or other HTS interfaces 
even if the material is clean in the bulk. There are at 
least two reasons for that. The first one is that quasipar- 
ticle reflection from realistic interfaces is diffusive, thus 
providing the isotropisation in the momentum space and 
the suppression of the d-wave component of the order 
parameter. The second one is the contamination of the 
material near interfaces as a result of fabrication process 
or electromigration in large scale application devices. As 
a result the formation of a thin disordered layer near HTS 
surfaces and interfaces is highly probable. An important 
question is whether or not superconducting correlations 
vanish in such a layer in the limit of small mean free path 
and what is the orbital structure of the superconducting 
state in the interface region. 

Surface peculiarities in d-w&ve superconductors were 
extensively discussed in the framework of the theoretical 
models based (Wii-specular quasiparticle reflection from 
clean interfaceaJrQ. Zero- and finite-bias anomalies pre- 
dicted in thesej-papers were recently observed experimen- 
tally in Refs.Ejta. In this paper we focus on the problem 
of the anomalous proximity effect between a c?-wave su- 
perconductor and a thin disordered layer in the limit of 



strong disorder. It is shown that an isotropic order pa- 
rameter is nucleated in such a layer even in the absence of 
the subdominant pairing interaction in the s-wave chan- 
nel. The spatially-resolved quasiparticle density of states 
is calculated. It is shown that zero- and finite-energy 
peaks are present in the surface density of states in the 
d-wave region. Zero-energy peaks are fully smeared out 
in the disordered layer, which is in a peculiar gapless su- 
perconducting state. 



II. PROXIMITY EFFECT AT THE INTERFACE 

Two approaches to the study of surface roughness ef- 
fects in unconventional superconductors were used pre- 
viously. In the first one it is assumed that the interface 
consists of facets with random orientations, compared to 
the crystallographic axes of the materials. According 
to the second approach, both sides of an ideal inter- 
face are coateci-by a so-called Ovchinnikov's thin disor- 
dered layerulljO. In the latter case the degree of disorder 
(or interface roughness) is measured by the ratio of the 
layer thickness d to the quasiparticle mean free path in 
the layer I. Up to now both approaches were used to 
study the smearing of Andreev surface bound states by 
weak disorder. Here we will concentrate on the regime of 
strong disorder. 

We consider the surface or a weakly transparent bar- 
rier in a d-wave superconductor oriented normal to the 
crystallographic ab plane. We assume that the crossover 
from the clean to the dirty limit takes place in a thin 
layer near the surface with mean free path £ and thick- 
ness d < y/^oZ, where £o is the coherence length of the 
bulk material. 

To study the proximity effect at the .witerface we use 
the quasiclassical Eilenberger equationsE3 with impurity 
scattering taken in the Born limit. For our purpose 
it is convenient to rewrite these equations in terms of 
functions $ + = (f{r,0) + f(r,6 + it))/ 2 and = 
f(r,0)-f(r,9 + 7T) 

4lu$ + + v cos 9^- = 4A 3 + - (g ($+) - $ + (g)) (1) 
dx T 

2vcos(9— ± = -(2u)+ - to))*. (2) 
dx T 

2 W cos0^ = (2A +-($+))$_ (3) 
dx T 
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w>0 



(X(9, 0')$-\ 



= 0. 



(4) 



Here g and / are respectively the normal and the anoma- 
lous quasiclassical propagators, u — irT(2n + 1) are the 
Matsubara frequencies, v is the Fermi velocity, x is the 
coordinate in the direction of the surface normal, 9 is the 
angle between the surface normal and quasiparticle tra- 
jectory, r = l/v and (...) = (l/27r) J^ n (...)d9. We assume 
that the Fermi surface has a cylindrical shape. 

For pure d-wave interaction the coupling constant may 
be written in the formtj 

X(9, 9') = X d (9, 9') = 2A cos(2(6» - a)) cos(2(0' - a)), 

where a is the misorientation angle between the crys- 
tallographic a axis and the surface normal. Then ac- 
cording to the self-consistency equation (^) angular and 
spatial dependencies of the pair potential are factorizcd 
A = \/2A(x) cos(2(6* — a)), i.e. A has pure d-wave angu- 
lar structure everywhere in the interface region. 

Far from the interface the bulk anomalous propagator 
also has the d-wave symmetry 



v^Aoo cos(2(0- a)) 
^ 2 + 2A^ cos 2 (2(9 -a))' 



(5) 



At the same time, as will be shown below, the angu- 
lar structure of the propagator & + (x,9) is essentially 
modified near the interface and an s-wave component of 
<&+(x, 9) is induced. 

To proceed further we have to supplement equations 
w hh the appropriate boundary conditions for the 
function and its derivative d<S> + (x)/dx at the in- 

terface between the clean and the disordered regions of a 
d-wave superconductor (at x = 0). These conditions can 
be derived by integration of the Eilenberger equations 
m a small region near the interface. In accor- 
dance with Refo, the first boundary condition is the 
continuity of $_ at the interface and can be written in 
the form 



£cos6» d$+(-0) ucos6»d$ + (+0) 



(<?(-0)) dx 



2^ 



dx 



(0) 



This condition manifests the current conservation across 
the interface. 

The second boundary condition depends on the 
backscattering properties of the interface. To account 
for such a backscattering we introduce a strongly disor- 
dered thin layer located near the interface at —8 < x < 0, 
which is characterized by the mean free path £$, where 
£g, 5 <C d,£. Assuming that all the boundaries are trans- 
parent, integrating Eq.(Q) in the interval —8 < x < 
and taking the limit 8 — > we arrive at 



dx 

where D = 28 /£ s . 



$+(+0)-4+(-0) 1 



(7) 



For D = Eq.(0) is a direct consequence of the 
continuity of the Eilenberger functions along quasiclas- 
sical trajectories which is valid for transparent SN- 
boundariesE-i. With increase of D the probability of 
quasiparticle penetration into the iV-layer decreases as 
D _1 , i.e. $+(-0) w L> _1 $ + (+0). It means that most 
quasiparticles are diffusively reflected back to the d-wave 
region at a length scale smaller than i. 

In the following we will consider the case-aLslrong dis- 
order £ <§; d. Then in accordance with Refs.Ej'lij it follows 
from the symmetry of the problem that the boundary 
condition at the totally reflecting free interface (x = —d) 
is 



d_ 

dx 



$+(-d) = 0. 



(8) 



Since £ <C d and d < \/£o£, the dirty limit condition 
£ <C £o is fulfilled in the disordered layer. It is straightfor- 
ward to show from Eqs.(Q), ([!]) that in this case the pair 
potential A in the disordered layer vanishes due to im- 
purity pair-breaking. Then the angle-averaged functions 
($+) and-4?) at — d < x <C — £ obey the dirty limit Usadel 
equationsLL3 in the form which formally coincides with the 
one valued for a normal metal with T cn = 0. Since the 
scale of variation of ($+) and (g) in this regime is of the 
order of the dirty limit coherence length VCo^, the func- 
tions (<&+) and (g) in a thin disordered layer d < y^ol 
are spatially-independent. 

As a result the Eilenberger equations in the region 
— d < x < are essentially simplified and have the solu- 
tion 



, cosh(fc(a; + d)) 

A — : — , $_ 

cosh(fcd) ' 



($+) cosh(fc(x + d)) 



= -2A 



sinh(fc(a: + d)) 
cosh(/cd) 

(9) 



<<?> 



cosh(fcd) 



($ + ) 2 + (.g) 2 = l. 



(10) 



where k = l/£ |cos 9\ . 

Making use of the boundary conditions Eqs.(^J), (Q) at 
x = and of Eqs.(^|), ([!(]), one can further reduce the 
problem to the solution of the Eilenberger equations (|l|)- 
(Q) in the clean d-wave superconductor (x > 0) 



dx 2 



dg_ 

dx 



A 
— < 



v |cos#| 
2uj ' 



to dx 



(11) 



(12) 



with the condition (j|) in the bulk (x 3> £o) an d the 
following boundary condition at x = 

{k (.9(0)) + D^} iU + (0) = $+(0) - ($+(0)) . (13) 
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FIG. 1. Angular dependencies of &+(x)/ttT c at different 
distances from the interface at T — 0.7T C . (a) Misorientation 
angle a — 0, (b) a = 5° 
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FIG. 3. Behavior of the pair potential near the interface 
for different misorientation angles a. Insert: dependence of 
{$+(0)} on a. 

In the vicinity of the interface (x < 0.3£o) the situation 
is just the opposite. In accordance with solution ([l4]), due 
to the angular dependence of k(9) oc |cos 6 1 ! the negative 
lobes are suppressed stronger than the positive ones, the 
function ($ + (a;)) changes sign to positive and reaches its 
maximum at x = 0. 

Note that for quasiparticle trajectories parallel to 
the surface, 9 = n/2, it follows from Eq.(|l4j) that 
$+(0,7r/2) = *(0,7r/2), while lim^ ±7r/2 $+(0, 6) = 
^(0, 7r/2) + X). This discontinuity is the manifestation of 
the simple fact that quasiparticles which propagate ex- 
actly parallel to the interface have information about the 
disordered region only via the local value of the pair po- 
tential, while for all other directions the direct interaction 
between both regions takes place. However, this discon- 
tinuity at 6 = ir/2 does not contribute to the result of 
the ang ular averaging of $ + in the boundary condition 
Eq.@. 

Fig. 3 shows the spatial variations of A(x) for different 
values of the angle a. As follows from Eq. (J14|) , the func- 
tion <I> + (0, 6*) near the interface has a contribution pro- 
portional to |cos 6*| cos(2(# — a)). This immediately leads 
to the result that the amplitude of the s-component in- 
duced into the disordered layer scales with misorientation 
angle a as ($+(0,q = 0)) cos 2a. At a — it/ 4 the super- 
conducting correlations are not induced into the disor- 
dered layer, i.e. ($+(0)) = 0. Further increase of a leads 
to a sigh change of the s-component. 

As is seen from Fig. 3, these qualitative considerations 
are in a good agreement with the results of exact numeri- 
cal calculations. In particular, for the d xy case (a = 7r/4) 



FIG. 4. The densities of states at T = 0.1T C in the disor- 
dered layer (solid triangles) and in the d-wave region at x — 
(open circles), x = £q (dashed), x = 2£o (dash-dotted) and in 
the bulk (dotted) 

it follows that ($+(0)) = 0. At the same time, it is 
worth mentioning that the pair potential at the inter- 
face, A(0,a = 7r/4), is nonzero, in contrast to the case 
of a specular reflecting boundary when A at a = 7r/4 
vanishes. In the case considered of diffusive surface scat- 
tering there is no symmetry requirement for the vanishing 
of $+(0,a = tt/4). 

In the whole temperature range the amplitude of the 
s-wave component ($+) induced into the disordered layer 
is an order of magnitude smaller compared to the ampli- 
tude of the order parameter in the bulk superconductor 
(see Fig. 2). That means that (3(0)) is close to unity for 
all temperatures. Thus, taking into account that (<?(0)) 
is independent of the Matsubara frequencies and that 
(0) ~ $+(0), we obtain from the boundary condi- 
tion Eq.fjll^) that at low temperature ($+(0)) oc u> for 
oj < A. As soon as u) exceeds the value of A, the function 
($+(0)) behaves as ($+(0)) oc lu~ 2 . The density of states 

N(e) = N(0)Re (g(0,e = iu)), where (g) = yfl - ($+) 2 
and N(0) is the normal state density of states. Therefore 
it follows from the property ( < I > +) oc w at small u> that 
the density of states N(e — Q)/Nq — 1, i.e. there is a 
gapless superconducting state in the disordered layer. 

To demonstrate this behavior explicitly we have cal- 
culated the density of states by numerical integration of 
Eqs.(ll)-(13) on the real energy axis making the substi- 
tution u> = —is in these equations. The results of calcula- 



tions of the normalized low-temperature density of states 
N(e)/N(0) = 1 at T = 0.1T C and a = are presented 
in Fig. 4. To take into account an inelastic scattering 
we have introduced the complex energy e = e + 17 with 
7 = 0.05T C . 

As is seen from Fig. 4, the density of states in the dis- 
ordered layer is gapless and has a rather weak and broad 
peak at an energy below the maximum value of the bulk 
pair potential A-bulk ■ This peak is a signature of the An- 
dreev bound states at finite energies. Note that there 
is no midgap (zero-energy) peak in the density of states 
in the disordered layer since there is no sign change of 
the order parameter in this region. At the same time 
the density of states at the surface of the d-wave region 
(at x = Q) exhibits a sharp zero-energy peak. In aceofc 
dance with the known results for specular interfacesnO'tJ 
this peak is due to the midgap states in the surface of 
the d-wave region . 

An important difference between specular and rough 
interfaces is that in the former case the midgap states 
occur only at nonzero values of the misorientation an- 
gle a 7^ 0, whereas in the latter case these states occur 
even for a = 0. This result has a simple physical inter- 
pretation. Due to the presence of the disordered surface 
layer incident and reflected quasiparticle trajectories are 
uncorrelated. Therefore for any a there is a finite prob- 
ability for an incident quasiparticle to experience a sign 
reversal of the pair potential upon reflection. Then the 
averaging over incoming trajectories yields the midgap 
states and as a result the nonvanishing zero-energy peak 
in the surface density of states appears. 

IV. CONCLUSIONS 

In conclusion, the proximity effect between a d-wave 
superconductor and a disordered surface layer is stud- 
ied theoretically in the regime of strong disorder. The 
boundary conditions for the Eilenberger equations are 
derived at the interface between the clean and the disor- 
dered regions. It is shown that superconducting correla- 
tions in the disordered layer do not vanish in the limit of 
small electronic mean free path and the isotropic super- 
conducting state is induced in such a layer. The crossover 
from this state to the d-wave pairing state in the bulk is 
studied by solving the Eilenberger equations. 

The quasiparticle density of states in the disordered 
layer is gapless and shows a broad peak at finite energies, 
while the density of states at the surface of the d-wave 
region exhibits a zero-bias peak due to sign reversal of 
the order parameter. This peak is fully smeared out in 
the disordered layer. 

The above phenomena have important consequences 
for the description of rough HTS interfaces and Joseph- 
son junctions, in particular for SNS junctions with a nor- 
mal metal interlayer. These effects will be discussed else- 
where. 
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